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Abstract 

It has recently been realized that supernova neutrino signals may be affected by shock propaga- 
tion over a time interval of a few seconds after bounce. In the standard three-neutrino oscillation 
scenario, such effects crucially depend on the neutrino level crossing probability Ph in the 1-3 sector. 
By using a simplified parametrization of the time-dependent supernova radial density profile, we 
explicitly show that simple analytical expressions for Ph accurately reproduce the phase-averaged 
results of numerical calculations in the relevant parameter space. Such expressions are then used 
to study the structure of Ph as a function of energy and time, with particular attention to cases 
involving multiple crossing along the shock profile. Illustrative applications are given in terms of 
positron spectra generated by supernova electron antineutrinos through inverse beta decay. 

PACS numbers: 14.60.Pq, 13.15.+g, 97.60.Bw 
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I. INTRODUCTION 



The realization that the supernova shock propagation can affect neutrino ffavor transitions 
Jl a few seconds after core bounce 01 is gaining increasing attention in the recent hterature 
on supernova neutrinos 0, B | H, 0,13, IE S Indeed, time-dependent variations of the 
neutrino potential in matter along the supernova shock profile can leave an interesting 
imprint on the energy and time structure of the (anti) neutrino signal at the Earth |8|, 
Tol |. Observation of such possible effects, although very challenging from an experimental 



viewpoint, would open a unique opportunity to study some aspects of the supernova shock 
dynamics in "real time." Moreover, the strong sensitivity of such effects to the neutrino 
oscillation parameters 0, 0, might, in principle, provide us with additional constraints 
(or hints, at least) on the neutrino mass spectrum and mixing angles. 

From a conservative viewpoint, it should be stressed that current simulations of core- 
collapse supernovae, despite remarkable efforts, may still require substantial physics im- 
provements l^l . It is not excluded that more refined simulations might significantly modify, 
e.g., the main features of shock profile discussed in ilj, and especially its propagation veloc- 



ity and density gradient, which govern the structure of the neutrino crossing probabilities 
0, H, Moreover, the time and energy dependence of the source (anti) neutrino fluxes 
might be more complicated and uncertain |l3l . l3 | than is customarily assumed, making it 
more difficult to identify shock-related signals. All these potential systematic uncertainties 
in the basic physics ingredients might be large enough to weaken the significance of shock- 
related signals, even with hypothetically large experimental statistics. However, despite all 
these caveats, the stakes are so high that further investigations on possible shock effects on 
supernova neutrinos are widely justified, in our opinion. 

The purpose of this work is to investigate in some detail the energy- and time-dependence 
of the neutrino crossing probabilities along the supernova shock. We assume the active 3iy 
oscillation scenario (Sec. II), and construct an empirical parameterization for the shock 
profile (Sec. III). Analytical calculations are shown to reproduce well (and to be much more 
convenient than) numerical approaches to the overall crossing probabilities, and are used 
to break down the multiple transition structure along the shock density profile, both in the 
energy domain (Sec. IV) and in the time domain (Sec. V). Finally, applications are given 
in terms of antineutrino event spectra observable at the Earth through inverse beta decay 
(Sec. VI). Conclusions and prospects for further work are given in Sec. VII. 

The bottom line of our work is that simple, analytical calculations of the crossing prob- 
abilities along the supernova shock profile appear to be adequate to study the structure of 
neutrino flavor transitions in the energy and time domain, for all practical purposes. The 
specific probability values estimated in this work, however, should not be taken too literally, 
being based on a simplified description of the matter density profile, which is subject to 
change as more detailed supernova shock simulations are performed and become publicly 
available. 



II. 3zv NOTATION AND ANALYTICAL APPROXIMATIONS 

Throughout this work, we consider only flavor transitions among the three known active 
neutrinos z/g, t'^, and i/^. In the following, we set the notation and briefly describe the 
analytical expressions for the crossing probabilities. 
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A. Notation for kinematics and dynamics 



The squared mass spectrum is parameterized as 

2"'+^'^^""') ' 

where the sign before Am^ distinguishes the cases of normal (+) and inverted (— ) hierarchy. 

The mixing matrix U is defined in terms of three rotation angles 6ij, ordered as for the quark 
mixing matrix jl5l |. 

U = U{9i2, 013,023). (2) 

When needed, we take 

6myeV^ = 7.3 x 10"^ , (3) 

sin^^ia = 0.315 , (4) 

corresponding to the so-called LMA-I best fit to the solar and reactor neutrino data jl6l |. 
We also take 

Am^ = 3 X 10^2 eV^ , (5) 

close to the best fit to the atmospheric and accelerator neutrino data detected in Super- 
Kamiokande 0|. The value of sin^ ^23 is basically irrelevant for supernova neutrinos, since 
it corresponds to an unobservable rotation in the (i/^, Ur) fiavor subspace. Concerning the 
mixing parameter sin^^ia, we will pick up several representative values, well below the 
current 3a upper limit (sin^ 6*13 < 0.05 17]). 

The mass gap hierarchy [dm? j /S.rn? <^ 1) and the smallness of sin^ 6'i3 guarantee, to a 
very good approximation, the factorization of the 3i/ neutrino dynamics into a "low" (L) 
and a "high" {H) 2v subsystem, whose dominant oscillation parameters are (5m^, ^12) and 
(Am^, 6*13), respectively (see |l8|, |l9| and references therein). The corresponding neutrino 
wavenumbers are defined as 

kL = 6m^/2E , (6) 
kn = Am^/2E , (7) 

where E is the neutrino energy. 

Potentially large matter effects jl^ are expected when either ki or ku take values close 
to the neutrino potential V in matter, 

V = V2GFNe{x) , (8) 

where A^e is the electron density at the supernova radius x. In appropriate units, 

-^ = 7.6xio-y;^, (9) 

eVVMeV g/cm^ ^ ' 

where Ye is the electron/nucleon number fraction (that we assume equal to 1/2), while p{x) 
is the radial mass density profile. We conventionally take fc^, ku, and V to be positive (as 
in Ref. ^^), irrespective of the neutrino type {v or v) and hierarchy (normal or inverted). 

The V{x) profile in supernovae is often approximated by a static power law, V{x) ~ a;"" 
with n ~ 3. In the presence of a shock wave, however, such static profile can be significantly 
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modified jH 0|- In particular the shock wave, while propagating outwards at supersonic 
speed, leaves behind a rarefaction zone, followed by a high-density region and by a sharp 
drop of density (down to the static value) at 

Xs = shock front radius . (10) 

This abrupt change of density occurs over a microscopic length scale governed by the mean 
free path of ions and electrons, which, for our purposes, can be effectively taken as zero.^ 
Correspondingly, a.t x = Xg the matter potential V{x) drops from 

V+ = \im_V{x) (11) 

to 

V- = lim^ V{x) , (12) 

with a typical ratio ^ 

e = — ^ 10 . (13) 

Given the peculiar step-like features of the matter potential at the shock front, in the 
following we discuss the analytical expressions of the crossing probabilities for x ^ Xg and 
X = Xg separately. We then study the general case of multiple crossings along the neutrino 
trajectory. 



B. Crossing probability at x ^ Xg 

As previously mentioned, potentially large matter effects are expected when kL,H ~ V- 
In particular, if the condition 

V{xh) = kn (14) 

occurs at some point xh{j^ Xg), there is a finite probability Ph of crossing between effective 
mass eigenstates in the H subsystem, which is given to a good approximation by the so-called 
double exponential formula |i2Q|] 



exp(27rr//A;HCOs2 6'i3) - 1 

= 77: 1—. (loj 

exp[27rrHkH) — 1 



(16) 



for both neutrinos and antineutrinos (see |19| and related bibliography therein), with 



d\nV{x) 



dx 



-1 



x=xh 

For all the V{x) profiles and sin^ ^13 values that will be used in the following, the above 
formula reduces to the well-known Landau-Zener (LZ) limit 

Ph ^ exp{-27TrHkH sin^ ^13) , (17) 



^ The authors of correctly recognized that the apparent smoothness of numerical shock front profiles 
0,1^ is actually an artifact of hydrodynamical simulations, which tend to widen the true (step-like) profile 
over several numerical grid zones. 
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at any v energy of practical interest. 

In principle, the occurrence of the condition fc^ = V{xi) at some radius xl could also 
induce neutrino level crossing in the L subsystem with probability Pl (generally different for 
neutrinos and antineutrinos; see, e.g., 0, H^l). However, for sin^^i2 as large as in Eq. (jl)) 
it turns out that the L-transition is adiabatic, 

Pi(z/) ~ ^ Pi(l7) , (18) 

at all points x ^ Xg. 



C. Crossing probability at x = 

The analytical approximations in Eqs. (fT7j) and ()18|) break down at x = x^, where the 
density gradient explodes (and th — > 0). In this case, the crossing probability Pg is de- 
termined by the conservation of flavor across the V{x) discontinuity, with the well-known 
result (see, e.g., [l8l|): 

P, = sin2(e-^^J , (19) 

where 9^ and ^~ are the effective mixing angles in matter immediately before (V = V~^) 
and after {V = V~) the shock front position Xg, respectively. For the H transition, such 
angles are defined by 

cos 20- ^ ^ cos 2^,3-1-^ A. 

'(cos 2^13 - V^/kn)^ + (sin 2^13)2 



for both neutrinos and antineutrinos.^ In the limit of very small ^13, the above equations im- 
ply Pg — 1 for kn G [V^^, V^], and P^ ~ otherwise. Notice that the extremely nonadiabatic 
limit Ps ~ 1 would also be obtained from Eq. (fT3j) for ^13 and th 0. However, the 
tempting approximation P, = 1 (for kn G [^~, V^]) does not work well for sin^ ^13 as large 
as O(10~^)-O(10~^), the "top-hat" behavior for Pg being significantly smeared out [i.e., a 
tail Ps 7^ is developed for kn outside the range [V~ , V~^], especially for kn ^ V~ , as easily 
understandable on the basis of Eq. ((201)] • Therefore, we adopt the exact Eqs. (fTI?|) and (p])) 
to calculate the H crossing probability at the shock front, without further approximations. 

Concerning the analogous crossing probability in the L sector, similar expressions apply, 
modulo the substitutions 6^13 — > 6*12 and k^ ki in Eqs. flT^ and ((201) The strong 
nonadiabatic character of the transition at the shock front can now lead to Pl 7^ in some 
corners of the parameter space [in contrast to Eq. (fTH|) ]. However, as we shall discuss at the 
end of Sec. IV C, such parameter values happen to have little phenomenological relevance. 
Therefore, unless otherwise noted, we shall formally take P^ ~ not only for x ^ Xg 
[Eq. (fTHj)]. but also for x = Xg. 



^ We remind the reader that, as proven in [l9|. Ph{v) = Ph{v), independently of the functional form of 
Ph- 

^ Contrary to the 7f-transition case, it is Pl{v) 7^ ^itEi general. Pl^P) can be obtained from Pl{v) 
through the further substitution /kh —V^/kh [iSj. 
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D. Analytical expressions for multiple crossings 



For Pl ~ 0, the formal expression for the Ve survival probability Pee |19l.l21| is exceedingly 
simple (up to Earth matter effects):^ 

{sin^ 6'i2 Pff (z/, normal), 
cos2 6'i2 (z7, normal), .^i) 
sin^ 6^12 (z/, inverted), 
cos^ $12 Ph 0^, inverted), 

where "normal" and "inverted" refer to the hierarchy type. In the above equation, we have 
neglected small additive terms of 0(sin^ ^13), which are not relevant for our discussion. From 
the above equations, it appears that Ph can modulate the (otherwise constant) survival 
probability of in normal hierarchy and of Ve in inverted hierarchy, thus providing an 
important handle to solve the current hierarchy ambiguity. 

In the presence of a propagating shock wave, the calculation of Ph is not entirely trivial. 
For the shock profiles studied in I3| , the condition in Eq. (|T^ can occur at up to three 
points Xs (i = 1, 2, 3): 

r(x,) = A;^^ (z = 1, 2, 3) , (22) 

In the most general case, two such points belong to the rarefaction zone [x\ < X2 < x^ 
and one to the static region below the shock front (xs > x^. The corresponding crossing 
probabilities are denoted as Pi, P2, and P3. The Pj's can be analitically calculated as in 
Eq. (fTTjl . with the inverse log-derivative [Eq. (fT6|l ] evaluated at the point Xj. In addition, 
one should consider the crossing probability P, at x = x^, as given by Eqs. (fT^ and ((201) • 
Therefore, the global Pu must be generally constructed in terms of four crossing probabilities 
Pi, P2, Ps, and P3, occurring at the points 

xi < X2 < Xs < X3 , (23) 

respectively. 

If we assume that all relative neutrino phases can be averaged away, and that the four 
transitions can be exactly factorized, the overall crossing probability Pu can be defined in 
terms of a simple matrix equation |18| 




=1.2, s, 3 



whose solution is 

P^^ = Pi + P2 + P, + P3 - 2(PiP2 + PiP, + P1P3 + P2P. + P2P3 + P.P3) 

+4(PiP2P, + P1P2P3 + P1P.P3 + P2P.P3) - 8PiP2P,P3 . (25) 

We will refer to this simple equation for our analytical calculations of Pu- 

Although the smallness of ^13 implies that each of the four transitions Pi,2,s,3 is (in 
general) well localized and can thus be separated from the others, the intrinsic accuracy of 
the factorization in Eq. (j24j) is not obvious a priori in the whole relevant parameter space. 



^ Examples of Earth matter effects will be illustrated in Sec. VI 
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There are zones of the shock profile where two transition points can become very close and 
eventually merge, or where the local density shape can be strongly different from the static, 
quasi-cubic power law which was used in to test the analytical recipe in of Eq. (fT3j) . 
Moreover, in the presence of multiple (interfering) transition amplitudes, it is worth checking 
that reasonable phase averaging effectively reproduces the incoherent factorization implicit 
in Eq. 

For such reasons, we think it is useful to compare the analytical calculations of Ph [based 
on Eqs. (fT^ and (fTTjl for x = Xi^2,3, and on Eqs. (fT^ and for x = Xg] with the results 
of a numerical (Runge-Kutta) evolution of the neutrino flavor propagation equations along 
representative shock density profiles. To our knowledge, such reassuring comparison has not 
been explicitly performed in the available literature.^ In order to do so in both the time and 
energy domain (Sees. IV and V, respectively), we need a parameterization of the shock wave 
profile which is (almost everywhere) continuous in both the radius x and in post-bounce 
time t. Unfortunately, a similar parameterization (or, equivalently, a dense numerical table) 
is not publicly available; only some representative profiles are graphically reported in Jl^ and 
j^.^ Moreover, only the profiles in exhibit (by construction) the correct step-like profile 
at x = Xg- In the following section, we thus introduce a simplified parametrization of the 
supernova density profile which is continuous in the {x, t) variables [excepting the shock- 
front discontinuity point]. Such empirical profile captures the main qualitative features of 
the time and radial dependence of the (supposed spherically symmetric) shock propagation, 
allowing us to perform a meaningful comparison of numerical and analytical results, and to 
break down explicitly the analytical calculation of Ph into its four components Pi, P2, Pg, 
and P3. Needless to say, such a simplified profile can and must be improved when more 
informative or refined supernova shock simulations become available. 

A final remark is in order. In this work we are mainly concerned with the factorization of 
transition within the H subsystem, since the factorization between the L and H subsystems is 
guaranteed for as low as in Eq. Q. However, small corrections to the L-H factorization 
might arise for values of Sm^ ~ few x 10^^ eV^, not considered here (see, e.g., for earlier 
discussions and jiol for more recent considerations). The analysis of these possible small 
corrections is beyond the scope of this work. 



III. EMPIRICAL PARAMETRIZATION OF THE SHOCK DENSITY PROFILE 

In this section we introduce a simplified, empirical parametrization of the shock density 
profile p, which reproduces the main features of the graphical profile in jl| and, at same 
time, is continuous in both x and t [excepting the case x{t) = Xg{t)]. In this way, numerical 
and analytical calculations of Ph can continuously cover the relevant parameter space. 

For post-bounce times t < 1 s, shock effects [l[ 1^ take place typically a.t V > kn and 
do not significantly affect the H subsystem. The density profile can thus be effectively 



^ The authors of Refs. and seem to have used a fully numerical and a mixed numerical+analytical ap- 
proach to the neutrino evolution equations, respectively. The authors of Ref. discuss only qualitatively 
the analytic form of P3, in the hypothesis tha,t_Pi ^ ^ P2- 

^ Notice that the bilogarithmic scales used in [l|, |3| make it difficult to recover potentially interesting shock 
details by graphical reduction. 



7 



approximated by its static limit po as taken from 



g/cm Vkm/ 



-2.4 



(26) 



For post-bounce times t ^ 1 s, the condition V ^ kn can be fulfilled, and the shock profile 
can thus modulate Ph- We characterize the profile in terms of the shock front position Xg 
and of its shape variation f{x) (with respect to po) for x < Xs- Formally, one can write 

t>ls p(x)=po(a:)-|f ^^"^^ ; x>S,' ^^^^ 
where ^ (the matter density enhancement a.t x = Xg) has been defined in Eq. ()13|) (see also 

Q)- 

In Eq. (j27j) . the function f{x) parametrizes the rarefaction zone ("hot bubble") above 
the shock front, characterized by a drop of density over more than one decade in x, and 
by an asymptotic density increase for x -C x^. After some trials, we have chosen the 
following (purely empirical) parametrization for f{x), which reproduces the main features 
of the graphical profiles in 

ln/(x) = [0.28 - 0.691n(x,/km)] ■ [arcsin(l - x/x,)]^'^ , (28) 



In the above equations we assume (following [l|) a slightly accelerating shock- front posi- 
tion Xs, 

Xs{t) =x^^+Vst+^ast\ (29) 
with parameters approximately given by 



X 







-4.6 X 10^ km , (30) 
Vs ^ 11.3 X 10^ km/s , (31) 
as ^ 0.2 X 10^ km/s^ . (32) 

Notice that: (i) The time dependence of p(x) in Eq. ^I7\ is imphcit, through the function 
Xs = Xs(t); and (m) at small radii, it is 

X < Xs p(x) ^ po(x) X 16 (xs/km)"^-^^ . (33) 

Figure ^ shows the neutrino potential V{x) as derived from the above parametrization 
[Eqs.(j2ni)-(|S21)] for representative post-bounce times > 1 s, as well as for the static profile 
(t ~ 0). These curves reasonably reproduce the main features of the shock profiles from 
[l|. The horizontal bands represent the ranges where the neutrino potential equals the two 
wavenumbers {kn = V and ki = V), for a representative energy interval E G [4, 70] MeV. 

As previously mentioned, a line at constant k^ can intersect the V{x) profile in Fig. ^ in 
at most three points x,, ordered as Xi < X2 (< x^) < X3, and corresponding to level crossing 
probabilities Pi, P2, and P3. A further crossing probability Ps is related to the step-like 
feature of the profile at the shock front x = x^. Two of the four critical points (xi, X2, x^, X3) 
can merge in the following cases: (1) at the bottom of the rarefaction zone (denoted as x^), 
where xi = X2; (2) at the shock front (x^), where it is X2 = Xs for = kn or X3 = Xs for 
V~ = kn- Single transitions can instead occur at early times or high energies {V{x) = kn at 
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X = only), as well as at late times or low energies {V{x) = kn aX x = xi only). Therefore, 
the three critical cases V{xr) = kn {xi = X2) and = ku (2:2,3 = Xg) are expected to mark 
significant changes in the behavior of Ph- 

We conclude this section with a few cautionary remarks. At present, the detailed shape 
of V{x) at the shock front is poorly known, since the physical requirement of a density 
discontinuity at implemented "by hand" as a necessary correction to the 

artificially smooth profile from simulations [8|.^ A side effect of the shock front steepening 
is the appearance of a local "cusp" at x = Xg (as apparent in Ref. [l| and even more in our 
Fig.Q), which might well be unphysical. However, we do not smooth it out, since it actually 
provides us with a useful, nontrivial check of numerical versus analytical calculations of Ph 
in an "extreme" condition (i.e., around a sudden change in the sign and value of the density 
gradient). Finally, we observe that some "leftover" turbulence might be expected in the 
rarefaction zone behind the shock front, inducing fuzzy variations of the local density. Such 
variations might lead to more than three solutions of Eq. (j22p in some cases, and thus to a 
more complicated and "random" structure for Ph. This possibility (not considered in this 
work) should be kept in mind when more refined supernova shock simulations will become 
available. 

IV. ANALYSIS OF CROSSING PROBABILITIES IN THE ENERGY DOMAIN 

In this Section we study in detail the energy dependence of Ph at fixed time. We pay 
particular attention to the discussion of multiple transitions along the shock density profile. 

A. Analytical vs numerical calculations of Ph{E) 

Figure 121 shows our calculation of Ph{E) at fixed post-bounce time t = 4s, for five 
representative values of sin^ ^13, ranging from 10^^ (top) to 10^^ (bottom). In the left panels, 
a direct comparison is made between analytical calculations (solid curves) and Runge-Kutta 
numerical calculations (dots), performed at 5 x 10^ equally-spaced points in the interval 
E G [0, 80] MeV. In general, up to four energy regimes can be identified, with significant 
changes at the three critical energies ~ 5, 50, and 67 MeV. These energies fulfill (for the 
V{x) profile at if: = 4 s) the three critical conditions kn = V'^, V~ , V{xr) which, as discussed 
in the previous section, signal the merging of two of the four possible level crossing points.^ 

In the left panel of Fig. |21 the results of the numerical calculations appear to be generally 
scattered around the analytical curve. The reason is that Runge-Kutta calculations keep 
track of the phase(s) of the level crossing amplitude(s), while our analytical approximations 
average out this information from the beginning. For a single numerical amplitude, the 
associated phase becomes ineffective after taking the squared modulus. For multiple (inter- 
fering) amplitudes, instead, the relative phases do provide an oscillating structure in Ph, 
which is apparent in the numerical results.^ 



^ In Ref. 1], the authors explicitly say to have steepened the shock front because "... in hydrodynamics 

calculations [it] may be softened by numerical techniques." 
^ Note that these critical conditions are independent of 6*13. 

^ In the left panels of Fig. |2 for the sake of graphical clearness, we show the spread of numerical results as 
scatter plots, rather than through rapidly oscillating curves. 
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In the right panels of Fig. |2l the oscillating structure of the numerical results is averaged 
out by a convolution with a Gaussian function (with one-sigma width of ±1 MeV), which 
simulates a generic "smearing" process (e.g., experimental energy resolution). The same 
convolution is applied to the analytical results. It can be seen that there is very good 
agreement between the two different (numerical and analytical) approaches after smearing. 
The results of Fig. |21 (and of other checks that we have performed for different post-bounce 
times) show that analytical calculations of Ph coincide with phase-averaged numerical results 
with very good accuracy, even close to the "critical" points where two transitions collapse. 
Moreover, analytical calculations are actually much more convenient than the numerical 
ones, which require a very high sampling rate in order to perform efficient phase-averaging 
and to prevent numerical artifacts. 

Let us now discuss in more detail the behavior of the analytical Ph in the left panels of 
Fig. 121 We remind that Pg [the crossing probability at the shock front, Eq. (jl9|) ] has almost 
a top-hat behavior: it rapidly drops from ~ 1 to zero for kn outside the range 
i.e., for E outside the range [5,50] MeV (for the specific profile used in Fig. |2)). This 
behavior depends only mildly on sin^ 6'i3 [as far as it is ^ O(10~^)]. The LZ probabilities 
-Pi,2,3 are instead exponentially suppressed as sin^ 6*13 increases [see Eq. [T3] . In particular, 
in the bottom panels of Fig. |21 (sin^ 6*13 = 10^^ and 10~^), Ph is dominated by Ps (with its 
characteristic top-hat behavior), with very little contributions from -Pi, 2,3- The dominance 
of a single crossing amplitude in Ph explains the suppression of the numerical oscillations 
in the same panels. 

Conversely, for sin^ ^^13 as small as 10~^ eV^ (top panel in Fig. |2I), the probabilities -Pi, 2,3 
(as well as Pg) can all become as large as ~ 1, rendering Ph ~ 1. In particular, for 
E ^ 5 MeV, only Pi ~ 1 is active (single crossing, well above the shock front). For 
5 ^ -E < 50 MeV, besides P, ^ 1 it is also Pi, 2 — 1 (two crossings in the rarefaction 
zone). For 50 ^ -E < 67 MeV, Pi, 2 are still active, P3 ~ 1 is switched on, and Ps drops to 
zero. Finally, for E ^ Q7 MeV, only P3 (level crossing well below the shock front) survives. 
In all such cases (characterized by either single or triple crossing with strong nonadiabatic 
character) one derives from Eq. ()25|) that Ph — 1- 

The intermediate panels in Fig. |21 corresponding to sin^ ^13 = 10^^ and 10~^, are not as 
easily understood as the previous ones, since at least one of the three probabilities Pi, 2,3 is 
definitely < 1. In such cases, Ph is better understood by explicitly separating its components 
P'i,2,s,35 as done in the next subsection. 

B. Ph{E) decomposition 

Figure 121 shows the decomposition of Ph{E) (analytical, unsmeared) into the four crossing 
probabilities Pi,2,s,3 for sin^ Oy^ = 10^^ (top panel) and sin^6'i3 = 10~^ (bottom panel). In 
both cases, the shock wave profile at t = 4 s (see Fig. (H) is assumed. The various transition 
probabilities are switched on and off at the three critical energies (-E ~ 5, 50, and 67 MeV) 
associated to such profile, as previously discussed. 

In the top panel of Fig. both Pi and P3 are strongly nonadiabatic (~ 1), and thus 
Ph = Pi ~ 1 and Ph = P3 ~ 1 in the single-transition regimes at low and high energy, 
respectively. At intermediate energies both P, and P2 are nonzero, but the nonadiabatic 
character of P2 is less pronounced, since at X2 (outer part of the rarefaction zone), the 
density gradient is smaller [and the scale factor in Eq. (jl6|) is larger] than at Xi or — 
at least in our simplified description of the density profile. Therefore, in the intermediate 
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regime 5 ^ E ^ 50 MeV, the strong nonadiabaticity of Pi and Ps implies Ph ~ P2 through 
Eq. (|25|) ■ In summary, in the top panel of Fig.|21 the overall crossing probability Ph initially 
takes the value Pi, then P2, and eventually P3, as the energy increases and passes through 
its critical values. All relevant parts of the density profile (static profile, shock front, and 
both sides of the rarefaction zone) can thus contribute to the Ph energy profile. 

In the bottom panel of Fig. El the higher value of sin^ 613 leads to an overall decrease of 
-Pi,2,3; particularly for the second transition (P2 ^ 0), while Pg remains strongly nonadiabatic. 
Therefore, the behavior of Ph is not obvious, except at very low energy (where only Pi > 0) 
or at very high energy (where only P3 > 0). We can only say that the smallness of P2 implies 
that the outer part of the rarefaction zone (just behind the shock front) are not relevant in 
this case. 

Summarizing, variations of sin^ ^13 for fixed density profile can modulate significantly the 
-Pi,2,3 contributions to the crossing probability Ph- The contribution of Pg remains instead 
stable and strongly nonadiabatic. As a consequence, different parts of the neutrino profile 
can play leading or subleading roles, according to the chosen value of ^13. When 6^13 will 
be known, one shall be able to gauge the relative importance of such contributions, and 
to determine which parts of the shock profile should be most intensively studied. In the 
meantime, all parts of the shock profile should be considered as potentially important and 
worth further study. 

C. Ph{E) at different times, and comments on 

We conclude the analysis in the energy domain by briefly discussing the variations of 
Ph{E) for different post-bounce times t. We also comment on the contribution of the shock- 
front discontinuity to P^. 

Figure |3 shows the function Ph{E) at three different times {t = 0, 4, and 8 s) for 
sin^ 6'i3 = 10~^ (top panel) and 10~^ (bottom panel). The curves at t = 4s have already 
been discussed. At earlier times t^ls(t = Osin Fig. the behavior of Ph becomes 
simpler, being basically dictated by a single LZ crossing along the static profile. At later 
times (t = 8 s), the critical energies at which kn = and kn = V{xr) move forward and 
partly fall outside the energy range relevant for supernova neutrinos. As a consequence, the 
nonadiabatic character of the innermost transition at Xi makes Ph ~ 1 in almost the first 
half of the energy range of Fig. El with a rapid drop to smaller values as the more adiabatic 
transition P2 becomes relevant. The drift of the zones where Ph is high or low has thus an 
interesting time structure, which is analyzed in more detail in the next section. 

We conclude this section with a few comments on the crossing probability Pl in the L 
subsystem. As mentioned in Sec. II C, the step-like nature of the shock front profile can 
make P^ ~ P^ 7^ for k^ roughly within the range [V^~,V^"'"]. In practice, however, this 
condition occurs only in the low-energy tail of the (anti)neutrino spectrum {E ^ few MeV) 
and at relatively late times (t > 10 s). Low-energy effects are suppressed by the lower cross 
section and by experimental detection thresholds, while late-time effects are intrinsically 
suppressed by the exponential decrease of the supernova neutrino luminosity on a timescale 
of a few seconds (see also Sec. VI). Moreover, the L-transition is never strongly nonadiabatic; 
in fact, from Eq. ()19p one can prove that Pl(z^) and PlO^) are always significantly smaller 
than cos^ 6^12 and sin^6'i2, respectively. For such reasons, effects related to P^ 7^ at the 
shock front are hardly observable, and have been neglected throughout this work. 
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V. ANALYSIS OF CROSSING PROBABILITIES IN THE TIME DOMAIN 



In this Section we study in detail the time dependence of Ph at fixed energy. As in 
Sec. IV, we start with a comparison between numerical and analytical calculations of Ph, 
and then we study the decomposition of Ph in terms of -Pi,2,s,3- 

A. Analytical vs numerical calculations of Pnit) 

Figure El shows our calculation of Pnit) at fixed neutrino energy E = 30 MeV and for five 
representative values of sin^ 613, ranging from 10~^ (top) to 10~^ (bottom). In the left panels, 
a direct comparison is made between analytical calculations (solid curves) and Runge-Kutta 
numerical calculations (dots), performed at 5 x 10'^ equally-spaced points in the shown time 
interval. Once again, the interference among multiple u transition amplitudes appears to 
generate fast oscillations, which scatter the numerical results around the analytical curves. 
In the right panels of Fig. 13 phase averaging (smearing) is instead enforced through a 
convolution with a "top-hat" time resolution function with ±0.5 s width. It can be seen 
that such phase averaging brings the numerical and analytical calculations in very good 
agreement at any t, up to small residual artifacts in the smeared numerical results, which 
would disappear by sampling the abscissa more densely (not shown). From the results of 
Figs. 121 and we conclude that simple, analytical calculations provide the correct phase- 
averaged value of Ph as a function of both E and t. Therefore, time-consuming numerical 
estimates of neutrino transitions along the shock profile 0, |^ can be accurately replaced by 
much faster and elementary calculations. This is one of the main results of our work. 

Let us now discuss in more detail the behavior of the analytical Ph in the left panels of 
Fig. In analogy with the previous discussion in the energy domain, the structure of Ph 
appears to be characterized by (up to) four different regimes also in the time domain. These 
regimes are separated by three "critical times," namely, t ~ 3.2 s (when kn equals V at the 
bottom of the rarefaction zone Xr), t ~ 3.3 s (when kn equals V~ at the bottom of the shock 
front), and t ~ 7.7s (when kn equals at the top of the shock front). [Of course, for 
E ^ 30 MeV such critical times will be different.] The probabilities -Pi, 2,5,3 play different 
roles in forming the global Pnit) in these different regimes. For t ^ 3.2 s, the relevant 
density profile is essentially static, and Ph — P3 oc const*^™ Similarly, for t ^ 7.7 s, it is 
Pff ~ Pi oc const^^'^ _ The single-transition values of Ph at early and late times are thus 
reduced by successive powers of ten as sin^ ^13 is reduced by factors of ten from top to bottom 
in Fig. El At intermediate times the situation is instead less obvious. For t ~ 3.2-3.3 s, the 
crossing probabilities P2 and Pg are switched on in succession, leading to a complex local 
structure in Pnit) (especially at small sin^ ^13); in this situation, neutrino matter effects are 
probing the lower part of both the rarefaction zone and of the shock front (see Fig. In 
the interval 3.3 ^ t ^ 7.7 s, Ps is always strongly nonadiabatic, P3 is negligible, while Pi, 2 
play a role only at small sin^ 6'i3 (being exponentially suppressed for increasing ^^13). For 
sin^ 6'i3 = 10~^ (bottom panels), only Pg is important in forming the global Ph- 

The most complicated cases in Fig. El thus emerge at "intermediate times" (i.e., when 
both the rarefaction zone and the shock front are being probed) and at "intermediate 6'i3 
mixing" (i.e., when Pi,2,3 are moderately nonadiabatic). In such decomposition of 

Ph into Pi, 2,5,3 is necessary to grasp the main features of Pnit)- 
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B. Pnit) decomposition 

Figure El shows the decomposition of -Pff(t) (analytical, unsmeared) into the four crossing 
probabilities -Pi,2,s,3 at = 30 MeV, for sin^ 6*13 = 10~^ (top panel) and sin^ ^13 = 10~^ 
(bottom panel). The various transition probabilities are switched on and off at the three 
critical times t ~ 3.2, 3.3, and 7.7 s. In both panels, P3 and Pi dominate at early and late 
times, respectively. At intermediate times, the structure of Ph depends sensitively [through 
Eq. (j2Sl)] upon the specific values of the crossing probabihties (-Pi,2) in the rarefaction zone 
{Ps being almost constant and close to 1). The results in this figure confirm that variations 
of can lead to dramatic variations in the role of each Pi (and thus of the different parts of 
the shock profile) in building the final crossing probability Ph as a function of both energy 
and time. 



VI. IMPACT ON OBSERVABLE POSITRON SPECTRA 

As an application of the analytical calculation of Ph described in the previous section, 
we study the effect of the shock propagation on the energy (and time) spectra of positrons 
detectable at the Earth through the inverse beta-decay reaction 

'T'e + P ^ n + . (34) 

We assume a "standard" supernova explosion at D = 10 kpc, releasing a total energy 
Eb = 3 X 10^^ erg, equally shared among the (anti)neutrino flavors u^, and distributed in 
time as 

L{t) = ^e''/^, (35) 
or 

where L is the luminosity for each z/^ flux, and the decay time is taken as r = 3 s 0| . 
For the sake of simplicity, we assume unpinched (normalized) Fermi-Dirac Ua spectra with 
time-independent temperatures T^, 

0.176 m 

47rD2 ei^/T. + 1 ' ^^^^ 

where we take = 4.5 MeV and = 6.5 MeV (x = /i, r). 1° 

In the presence of oscillations, the supernova z/g spectrum is given by 

f{E) = h{E)P,, + h{E){l - Pee), (37) 

where Pee is taken from Eq. (pT|l (antineutrino case). The convolution of f{E) with the 
differential cross section da{E, E^o^) I dE^o^ il^ provides the event spectrum dN/dE^os in 
terms of the positron energy Epos- We assume perfect energy resolution and zero threshold, 
in order to show the shock effects in the most favorable conditions. It is understood that 
such effects will be somewhat degraded in realistic cases, depending on detector details. 
Finally, just to fix the overall scale, we assume a detector volume corresponding to 32 kton 
of water. However, for the main purposes of our discussion, the event rate units could be 
taken as arbitrary in the following figures. 

10 Pqj, recent discussions of supernova neutrino energy spectra, see 
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Figure [7| shows our calculated positron spectra for no oscillations (Pee = 1, dashed curve) 
and hypothetical full conversion (Pee = 0, dot-dashed curve). For partial conversion, Eq. ()37|) 
implies that the final positron spectrum is basically a linear combination of the previous 
two spectra. Examples are given in Fig. [7|for sin^ ^13 = 10~^, for both normal hierarchy 
(dotted curve) and inverted hierarchy (solid curve). For normal hierarchy, the coefficients 
of the combination are energy- independent, namely, Pge = cos^^i2 and 1 — Pee = sin^ ^12 
[see Eq. |^ . Conversely, for inverse hierarchy, such coefficients become energy- and time- 
dependent through Ph [Pee = cos^ 9i2Ph{E)]. In the specific example of Fig. [71 the time 
is fixed through integration over a representative bin [t = 6 ± 0.25 s]. The energy profile 
of Ph{E) then governs the relative weight of the spectra fe and fx in the combination. As 
Ph decreases from its low-energy value at the first critical energy (see Figs. El and 01 and 
related comments), the weight of increases and the positron spectrum is shifted to higher 
energies. The "shoulder" around the peak of the solid curve is thus an imprint of the shock 
passage. This shoulder is expected to drift for time bins different from t = 6 ± 0.25 s. 

Figure IHl shows a set of successive positron energy spectra, binned in time intervals of 0.5 s, 
from about 6 to 9 seconds after bounce. The overall decrease of the spectra for increasing 
t is due to the decrease in luminosity. The top panel does not include Earth effects, while 
the bottom panels includes representative Earth mantle crossing effects, calculated in the 
same conditions as in As in Fig. [71 the cases of normal and inverted hierarchy are 

distinguished through dotted and solid curves, respectively, and sin^ ^13 is fixed at 10~^. The 
top panel shows how the shock passage can leave an observable, time-dependent spectral 
deformation, in the case of inverse hierarchy. This possibility is very exciting, although one 
might realistically hope to see no more than variations in the first spectral moments j^, E3| ■ 
From the comparison of the top and bottom panel in Fig. IHl it turns out that the spectral 
modulation due to the shock might be partly hindered, at late times, by the additional mod- 
ulation generated by Earth effects. In other words, the "wiggles" produced by oscillations 
in the Earth might make it difficult to identify the position of the genuine, shock- dependent 
spectral deformations at different times. Therefore, in some cases Earth effects may not 
necessarily represent an additional handle fld\ to study the shock dynamics via neutrinos. 




Finally, Fig. 1^1 shows the drop in the positron event rate (arbitrary units) as a function 
of time and for four representative values of sin^ 613, ranging from 10~^ to 10^^. Four rates 
are displayed in each panel, integrated over consecutive 10 MeV energy bins. For normal 
hierarchy (dotted curves), the rates decrease according to the assumed exponential law in 
Eq. ()35p. For inverted hierarchy (solid curves), the characteristic time structure of Pnit) 
(see Fig. is instead reflected in a characteristic deviation of the rate decrease with respect 
to a pure exponential drop: The lower Pnit), the higher the rate (enriched in events) for 
inverse hierarchy, as compared with normal hierarchy. The two hierarchies become instead 
indistinguishable when Pff(t) ~ 1 (see Fig. EJ, since Pee — cos^ 6'i2 for both normal and 
inverted neutrino mass spectra [Eq. (pT|l ]. Roughly speaking, the time structure of the event 
rates in Fig. El reflects (although "up-side-down") the Pff(t) pattern in Fi^. El 

Shock- wave signatures in the time domain have been emphasized in [ij. It goes without 
saying that very high statistics would be needed to discriminate such time structures in 
real experiments.^^ Nevertheless, the reward could be very high. Let us consider, e.g., the 



A more complete and quantitative study of the observability of shock-induced structures in supernova 
(anti)neutrino signals will be performed elsewhere, in the context of future large-volume, high-statistics 
detectors. 
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"bathtub" pattern of the sohd curves in the lower right panel of Fig. El (sin^6'i3 = 10~^). 
The (hypothetical) experimental detection of such rather distinctive pattern would, at the 
same time: (i) provide a snapshot (if not a "movie") of the shock wave propagation; (ii) 
prove that the neutrino mass hierarchy is inverted; and [iii) put a significant lower bound 
on sin^^i3. Each of these results would have a dramatic impact on our understanding of 
both supernova and neutrino properties. 



VII. SUMMARY AND PROSPECTS 



Supernova shock propagation can produce observable effects in the energy and time struc- 
ture of the neutrino signal, as suggested in jl[ and also investigated in 0, In the cal- 
culation of such effects for 3z/ oscillations, the neutrino crossing probability Ph in the 1-3 
neutrino subsystem plays an important role. The evaluation of Ph can be done either nu- 
merically or analytically, for given shock profiles in space and time. By assuming a simplified 
description of the time- dependent density profile, we have explicitly shown that simple an- 
alytical calculations accurately reproduce the phase-averaged numerical results, also in the 
nonobvious case of multiple transitions along the shock profile. The analytical approach 
has then been used to study some relevant characteristics of the structure of Ph in energy 
and time. This structure can in part be reflected in observable signals, as we have shown 
through some simple and selected examples. 

The simplifications used in this work, and especially the empirical parameterization of 
the density profile, do not alter our main conclusions about the validity and usefulness of 
a simple analytical approach to Ph, as compared with brute- force numerical calculations. 
However, such simplifications may affect the specific values of Ph that we have estimated at 
given energy and time. These values (and the corresponding observable effects) should not 
be taken too literally, since they depend on features of the shock profile which are beyond 
our control. In particular, the shape of the rarefaction zone left behind by the shock wave 
should be constrained, if possible, by dedicated supernova simulations with higher resolution 
in space and time. When such issues will be clarified, more realistic calculations of Ph, and 
of its imprint on observable neutrino signals at the Earth, will become possible. 
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neutrino potential profile 
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FIG. 1: Neutrino potential y as a function of the supernova radius x, for different values of 
the post-bounce time t, as derived from our simplified parametrization of the profiles in P|- The 
horizontal bands represent typical ranges of the neutrino wavenumbers kn and ki,. See the text 
for details. 
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FIG. 2: Comparison of numerical and analytical calculations of the crossing probability Ph{E) 
(dots and solid curves, respectively) at t = 4s and for five representative values of sin^ ^13. In the 
right panels, phase effects are energy- averaged. 
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Ph components (t = 4.0 s) 
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FIG. 3: Ph{E) in terms of its components -Pi,2,s,3 at t = 4s, for two representative values of 
sin^ (9i3. 
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Ph(E) at different times 
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FIG. 4: Ph{E) at t = 0, 4, and 8 s, for two representative values of sin^ 6 
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FIG. 5: Comparison of numerical and analytical calculations of the crossing probability Pnit) 
(dots and solid curves, respectively) at -E = 30 MeV and for five representative values of sin^^ia. 
In the right panels, phase effects are time-averaged. 
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Ph components 
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Energy spectra for u^(p, n) e"^ for t = 6.0±0.25 s 
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FIG. 7: Representative positron energy spectra of supernova antineutrino events from inverse 
beta decay, for different degrees of Vg conversion into F^- Daslied curve: No conversion. Dot- 
dashed curve: Full conversion. Dotted line: Partial conversion for normal hierarchy. Solid line: 
Partial conversion for inverse hierarchy (calculated for sin^ ^13 = 10~^, and integrated over the time 
interval t = 6 it 0.25 s). The latter case shows the shock imprint as a "shoulder" in the spectrum 
range E ~ 10-20 MeV, . 
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Energy spectrum for u^(p,n) e'^ for different times 
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FIG. 8: Positron energy spectra for successive time bins, calculated for sin^ ^13 = 10~^ and for 
normal and inverted hierarchy (dotted and solid curves, respectively). In the case of inverted 
hierarchy, shock-induced spectral deformation move forward in energy for increasing post-bounce 
time. The bottom panel includes representative effects of Earth mantle crossing, which induce 
additional spectral deformations. 
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FIG. 9: Time dependence of the positron event rate (binned in specified energy intervals) for 
four representative values of sin^^is. For normal hierarchy (dotted lines), the event rate is just 
proportional to the neutrino luminosity (assumed to decay exponentially in time). For inverted 
hierarchy (solid lines) an additional time structure is induced by the shock propagation. 
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